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The aim of this paper is to employ some estimation methods, namely, modified moment method and standard 
Bayes method to estimate the two parameters of the first kind Pareto distribution and comparing the performance of these 
two methods by using the Mont Carlo simulation technique. Moreover, some Applications related to Pareto distribution, in 
particular, human population statistics and tropical cyclones in Japan are given. The Chi Square goodness of fit test is 
presented for each case. 
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1. INTRODUCTION 

In this paper, the first kind Pareto distribution Par(a,c) is studied Pareto distribution is widely used as a model in 
many areas of application, for instance, the values of oil reserves in oil fields, sizes of sand particles, areas burnt in forest 
fires, economics, human population statistics, rainfall, and others. Professor of economics Vilfredo Pareto (1848-1923) 
was the first who proposed this distribution when he tried to put a comprehensive law to deal with the distribution of 
income for a particular community[l]. The remaining of this paper is organized as follows: in section 2, basic concepts on 
Pareto distribution are given, in section 3, two estimation methods are presented. Namely, the modified moment method 
and standard Bayes method. Monte Carlo results including the estimated values of the parameters as well as the estimation 
of some properties such as bias, mean, variance, kurtosis, skewness and MSE are presented in section 4. In section 5 we 
state the Chi Square goodness of fit test. Some real life applications related to Pareto distribution are presented in section 6. 
In section 7 the discussion of the results and conclusions are given. 

2. BASIC CONCEPTS ON PARETO DISTRIBUTION:^] 

A continuous random variable X is said to have Pareto distribution, denoted by X~Par(a,c) if it has Probability 
density function is: 


ABSTRACT 


f(x;a,c) = s x «+i 

0 


x>c ,c> 0, a>0 


( 1 ) 


otherwise 


Where a ,c are the shape and scale parameters respectively. 
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The greater value can be obtained for the function above, if x=c: 



C 


Pareto distribution depends on two parameters a and c, a wide variety of distribution shapes and seals can be 
generated from suitable choices of a and c, figures (1) show some Pareto p.d.f’s: 



The cumulative distribution function (cdf) ofX is: 

F(x;a,c) = P{X < x) 

X 

= |/(m)c/m 


F(x;a,c) = 1 - 



( 2 ) 


3. ESTIMATION METHODS 

3-1 The Modified Moments Method: [3] 

Let X l ,X 2 ,...,X n be a random sample from par (a, c) and let y! be the first order statistic of the sample, then 
the p.d.f of yds given as: 



Uccc na y; (na+1) , 


c — Ti 00 

e.w. 


Hence: 
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E b\) 


nac 
{na- 1) 


Equating y x to E(y,) in above: 


nac 

y 1 = — — r 

na -1 


Moreover, it can be shown that: 



a - 1 


Solving the last two equations yields the following estimators: 

x _ nX - y j 


na 


- y 1 ) 
~ _ {na- l)>’i 


(3) 

(4) 


3-2 The Standard Bayes Estimator: [4] [5] [6] 

Bayes style in estimation depends on the assumption that the Prior information available on the parameter to be 
estimated, can be formulated as the probability density function which is called the Prior probability density function (Prior 
P.d.f), and Bayes estimator depends on the Posterior probability density function (Posterior. P.d.f): 

h(e I x, x n ) oc L(x x x n ; 0)J(O) 


Where, L(x x ,..., X n ',0) is the Likelihood function 
J(6) - Prior P.d.f for 0 

h(e I Xj x n ) = kL(x l ,..., x n ; 0)J(0) 

k~ l = | L(x l ,...,x n ;0)J(0)d0 ; k is cons tan t 

V0 

0g =E(0/X l ,...,Xj 

For the case of Pareto Distribution and by using non informative Prior (p.d.f): 
According to Jeffrey method we have: 

JAa) oc — ,a>0 

1 a s 
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J 2 (c) oc — , c > 0 

c' 

J(a, c) oc J x (a) • J 2 (c) 

Where, each of 7, (a) ,J 2 (c) is the prior p.d.f for a,c respective ly 

J(a,c)cc -^-y 
a c 

Let Xi, X 2 ...,x„ be a random sample from Par(a,c), then the likelihood function is: 

n 

L(x 1 ,x 2 x n ;a,c) = n f {x t l a,c) 

i = 1 

" ac a 

L(x 1 ,x 2 ,.., 9 x n ;a,c ) = n 

i=I x,. 

we get: 


h(a,d x,, 

h(a,c/ x lv 
Let r=l 
h(a,c/ x, , 


-,W) 0C 


1 a"c” 


n/ ' aV Na+1 

IT X; 

7=1 


,x„) cc a"- s c" a ~ r exp 


-(« + l)X lnx / 


..,x n ) = k a n - s c na ~ l exp 


- (a + l)X lnx - 


k~ x =\ J a n ~ s c na ~ x exp 


-(a + Dl'n-v, 


|7c7a 



( " 


na 

71(11 

f 7 , X 

exp 

-E ln x, 

) 

C 

a ns exp 

-«Elnx,- 


V i=l 

A o 

na 

0 

V 1=1 ) 


= exp 


C „ X"° Tiff / 

-Ylnx, fa-^exp 
V t=i Ji na 


-a^lnx 

V 7=1 ' J 


da 


exp -^ku 


oo 

( n V 

exp 

-a\ -n lnx (1) 

0 

V i=i )_ 


da 
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exp 


V 1=1 


oo 

n 

r M 

f n-^-l 

1 a exp 

-a]Tln 

x t 

0 

1=1 

l x (p JJ 


k~ X 


r(„-o 


" 

( w 

E ln 

Xj 

/= i 

l x (l) J_ 


l>(f' 

n |_ »=i V x d) y 

f Vi 1 T(n-S) 

eX P _ ™ x i 


na" V“ 'exp -oT lnx ( 


h(a,c\ X l ,...,X n ) - 

\ t=l 

r „ f \ 

n-(«- S ) 

r (n-s) 

Zb — 



{ x mJ 


c = x (l) = min{X i ,X 2 ,. 



M(a X l ,...,X n ) = ^h(a,c 

1 X u ...,X n )dx 


C 


(5) 


M(a\X l ,...,X„)= 


a n 

-«E ln W 

_1 e - 


)(n-s) 


n—s 


w V c ) 



Suppose the Loss Function is Squared loss function then Standard estimator for shape parameter OC is mean of 
posterior distribution: 


r 


a = E 


a 

\x 



If s=0 then the standard Bayes estimator (X is the MLE, and if s=2 then Bayes estimator (X is minimum variance 
unbiased estimator (MVUE) of a 


a = 


n —2 


5> 




\C J 


( 6 ) 


4. MONTE CARLO RESULTS 

In this section, a large scale of samples generated from the Pareto distribution (par(3,l)) with size n(m); n=size of 
sample and m=Iterations, 5(1000), 10(1000), 20(1000), 30(1000), 50(1000), 75(1000), 100(1000), 250(1000), 500(1000) 
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and 1000(1000), These samples are used to estimate [OC,C Jof the Pareto distribution by using two method of estimator, 
Modified moments method and Standard Bayes method, and estimation some properties such us bias, mean, variance, 
kurtosis, skewness and MSE. 

In table (1), used equation (3), (4), (5) and (6) by using Matlab R2013b we get a simulation of the estimator 
(i a,c ), are displayed in table (1): 


Table 1: Parameters Estimation 


Sample size n( m) 

Estimation of (a,c) 

M.M.M 

B.M 

5(1000) 

(4.1729, 1.0077) 

(2.9343, 1.0713) 

10(1000) 

(3.4438, 1.0015) 

(2.9240, 1.0346) 

20(1000) 

(3.2797, 1.0000) 

(3.0196, 1.0164) 

30(1000) 

(3.1681,0.9995) 

(3.0063, 1.0106) 

50(1000) 

(3.1150, 1.0004) 

(3.0079, 1.0070) 

75(1000) 

(3.0870, 0.9997) 

(3.0128, 1.0041) 

100(1000) 

(3.0540, 0.9999) 

(2.9944, 1.0032) 

250(1000) 

(3.0341, 1.0001) 

(3.0094, 1.0014) 

500(1000) 

(3.0057, 1.0000) 

(2.9965, 1.0007) 

1000(1000) 

(3.0065, 1.0000) 

(3.0007, 1.0003) 


Table(l) show that methods M.M.M and B.M give a good agreement between the true values a = 3 and C = 1 
with the estimators (a) and (c) for all sample sizes, but the M.M.M give small values for c and B.M the better To 
estimation (X for all sample sizes. The following are some Pareto figures where par(3,l) for different sample sizes: 
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The biases of estimators CC and C for Pareto distribution which can be found by: 


bias (a) = E{a)~ a 
bias(c) = J £’(c)-c 

are displayed in tables (2) and (3): 


Table (2): Bias of Estimator (a) 


Sample size n( m) 

Bias of Estimation (a ) 

M.M.M 

B.M 

5(1000) 

1.1729 

-0.0657 

10(1000) 

0.4438 

-0.0760 

20(1000) 

0.2797 

0.0196 

30(1000) 

0.1681 

0.0063 

50(1000) 

0.1150 

0.0079 
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75(1000) 

0.0870 

0.0128 

100(1000) 

0.0540 

-0.0056 

250(1000) 

0.0341 

0.0094 

500(1000) 

0.0057 

-0.0035 

1000(1000) 

0.0065 

6.6286e-04 


Table 3: Bias of Estimator (c) 


Sample size n( in) 

Bias of Estimation (c) 

M.M.M 

B.M 

5(1000) 

0.0077 

0.0713 

10(1000) 

0.0015 

0.0346 

20(1000) 

4.3803e-05 

0.0164 

30(1000) 

-4.7243e-04 

0.0106 

50(1000) 

4.3246e-04 

0.0070 

75(1000) 

-2.6800e-04 

0.0041 

100(1000) 

-7.5317e-05 

0.0032 

250(1000) 

7.3703e-05 

0.0014 

500(1000) 

2.81 17e-05 

6.9569e-04 

1000(1000) 

-9.2540e-06 

3.2393e-04 


Table (4): show the variance of estimator (a ) which can be obtained by using the equation: 

2 

Var(a ) = — ^ («, - « ) 

n i=i 


Table 4: Variance of Estimator (d) 


Sample size n(m) 

Variance of Estimation (d ) 

M.M.M 

B.M 

5(1000) 

6.3794 

3.6320 

10(1000) 

1.3521 

1.0715 

20(1000) 

0.5802 

0.5058 

30(1000) 

0.4114 

0.3687 

50(1000) 

0.2199 

0.1938 

75(1000) 

0.1579 

0.1399 

100(1000) 

0.1023 

0.0921 

250(1000) 

0.0412 

0.0345 

500(1000) 

0.0227 

0.0173 

1000(1000) 

0.0123 

0.0098 


Table (4) show that the variance values of the estimator (a) by using the B.M are smaller than the variance values 


of M.M.M. 
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Table (5): shows the variance of estimator (c) which can be obtained by using the equation: 

Kar(c) = -J(c,-c) 

n ;=i 


Table 5: Variance of Estimator (c) 


Sample size n 

Variance of Estimation (c) 

M.M.M 

B.M 

5(1000) 

0.0054 

0.0052 

10(1000) 

0.0013 

0.0012 

20(1000) 

2.7613e-04 

2.6987e-04 

30(1000) 

1.2335e-04 

1.2494e-04 

50(1000) 

5.1961e-05 

5.1353e-05 

75(1000) 

1.8262e-05 

1.8186e-05 

100(1000) 

9.9469e-06 

9.8536e-06 

250(1000) 

1.9483e-06 

1.9391e-06 

500(1000) 

4.8348e-07 

4.8337e-07 

1000(1000) 

1.0994e-07 

1.0966e-07 


Table (5): show that the variance values of the estimator (c) by using the M.M.M are greater than the variance 
values of B.M. 


The m.s.e of estimator [OC ) and ) can be obtained by: 

ms.e(a) = Vat{a\ + [fo'a^a)] 2 
m.s.e(c) = Vai{c) + [/?/as(c)] 2 

Tables (6) and (7) show m.s.e of (a) and (c) by the two estimations methods. 


Table 6: Mean Square Error of Estimator (a) 


Sample size n(m) 

Mean Square error of Estimator (a) 

M.M.M 

B.M 

5(1000) 

7.7551 

3.6363 

10(1000) 

1.5490 

1.0773 

20(1000) 

0.6584 

0.5062 

30(1000) 

0.4397 

0.3688 

50(1000) 

0.2332 

0.1938 

75(1000) 

0.1655 

0.1400 

100(1000) 

0.1052 

0.0921 

250(1000) 

0.0423 

0.0346 

500(1000) 

0.0227 

0.0174 

1000(1000) 

0.0124 

0.0098 
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For all sample sizes in table (6) the m.s,e of cc given by B.M is less than other method results. 


Table 7: Mean Square Error of Estimator (c 


Sample size n(m) 

Mean Square Error of Estimator (c) 

M.M.M 

B.M 

5(1000) 

0.0054 

0.0103 

10(1000) 

0.0013 

0.0024 

20(1000) 

2.7613e-04 

5.3748e-04 

30(1000) 

1.2357e-04 

2.3677e-04 

50(1000) 

5.2148e-05 

1.0097e-04 

75(1000) 

1.8334e-05 

3.5333e-05 

100(1000) 

9.9525e-06 

2.0394e-05 

250(1000) 

1.9538e-06 

3.8988e-06 

500(1000) 

4.8427e-07 

9.6735e-07 

1000(1000) 

1.1003e-07 

2.1458e-07 


For all sample sizes in table (7) the m.s.e of C given by M.M.M is Less than B.M results. 
The skewness of estimators (a) and (c) which can be obtained by: 


Sekwensfqa) 


1 

n 


£(a,.) 3 -3a^(a ; ) 2 +3 <f £(«, )- a 3 



Sekwensdf) 


\_ 

n 





Are shown in tables (8) and (9): 


Table 8: Skewness of Estimator (a ) 


Sample size n(m) 

Skewness of Estimator (a) 

M.M.M 

B.M 

5(1000) 

1.8293e+03 

1.6368e+03 

10(1000) 

2.7416e+03 

2.4061e+03 

20(1000) 

4.0300e+03 

3.8694e+03 

30(1000) 

4.0474e+03 

4.0784e+03 

50(1000) 

5.8659e+03 

6.3890e+03 

75(1000) 

6.2495e+03 

6.9717e+03 

100(1000) 

8.7039e+03 

9.5994e+03 

250(1000) 

1.3363e+04 

1.7009e+04 

500(1000) 

1.5914e+04 

2.3534e+04 

1000(1000) 

1.9849e+04 

2.8032e+04 
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Table 9: Skewness of Estimator (c) 


Sample size n(m) 

Skewness of Estimator (c) 

M.M.M 

B.M 

5(1000) 

5.1883e+05 

6.4987e+05 

10(1000) 

2.2238e+06 

2.6022e+06 

20(1000) 

1.0887e+07 

1.1829e+07 

30(1000) 

2.4273e+07 

2.4610e+07 

50(1000) 

5.3412e+07 

5.5448e+07 

75(1000) 

1.7054e+08 

1.7389e+08 

100(1000) 

3.1837e+08 

3.2614e+08 

250(1000) 

1.4861e+09 

1.4861e+09 

500(1000) 

5.9439e+09 

5.9577e+09 

1000(1000) 

2.7404e+10 

2.7538e+10 


Tables (8) and (9) show that the two methods give a very small skewness to left and to the right which indicate 
that the estimators (<£ ) and (f ) approach rapidly to normal distribution. 

The kurtosis of estimators (a) and (c) which can be obtained by: 


kurtosis(a)= — 


Z(“< ) 4 - ) 3 + 6 « 2 2>, )- 4 « 3 Z(«. ) 

i=l i = 1 i=l i=l 

TT 


-3 


kurtosis(c) = — 


i=i ;=i ;=i /=i 

V 


-3 


Are shown in Tables (10) and (1 1): 


Table 10: Kurtosis of Estimator ((5) 


Sample size n(m) 

Kurtosis of Estimator (<f ) 

M.M.M 

B.M 

5(1000) 

-1.6462e+03 

1.1048e+03 

10(1000) 

-4.5059e+04 

-3.5588e+04 

20(1000) 

-9.4167e+04 

-8.6615e+04 

30(1000) 

-1.0597e+05 

-1.0489e+05 

50(1000) 

-2.0270e+05 

-2.2380e+05 

75(1000) 

-2.5027e+05 

-2.8645e+05 

100(1000) 

-4.2438e+05 

-4.7937e+05 

250(1000) 

-1.0100e+06 

-1.3876e+06 

500(1000) 

-1.5968e+06 

-2.6849e+06 

1000(1000) 

-2.6974e+06 

-4.2693e+06 
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Table 11: Kurtosis of Estimator (c ) 


Sample size n(m) 

Kurtosis of Estimator (c ) 

M.M.M 

B.M 

5(1000) 

-7.6843e+06 

-1.3735e+07 

10(1000) 

-6.3665e+07 

-9.323 le+07 

20(1000) 

-6.5712e+08 

-8.0445e+08 

30(1000) 

-2.1821e+09 

-2.3687e+09 

50(1000) 

-7.4419e+09 

-8.1295e+09 

75(1000) 

-3.9878e+10 

-4.2003e+10 

100(1000) 

-1.0100e+ll 

-1.0637e+l 1 

250(1000) 

-1.0545e+12 

-1.0787e+12 

500(1000) 

-8.5586e+12 

-8.6195e+12 

1000(1000) 

-8.2726e+13 

-8.3434e+13 


5. CHI-SQUARE GOODNESS-OF-FIT TEST:[5| 

Pearson Chi-Square goodness of fit test can be used for testing the null hypothesis involving discrete as well as 
continuous distribution, the Pearson chi-square test uses the density function of the random variable X. 

Let Xj , X 2 , . . . , X n be a random sample from a random variable X with probability density function f(x). We wish 
to test the null hypothesis 

H 0 -X~ f(x) 

Against: 

Hf.Xf f (x) 

For this we gusset the n observations into k mutually exclusive cells and we let p t be the probability that one x.of 
the sample is fall in cell i, Y j be the observed number of the cell i, and np i be the expected number of the cell i, i=l,2,...,n. 
Since the sample size is n, the number of observations expected to fall in the ith interval is equal to tip f . Then: 

8 =^ 

r-A tip i 

has an approximate chi-square distribution with k -1 

We reject the null hypotheses Ho if Q > X\- a (A 
Square distribution with l — a level of sig. and (k-1) degrees 

6. SOME REAL LIFE APPLICATIONS RELATED 
6-1 Application (1): Human Population Statistics 
6-1-1 Sweden 

In this application we use the population statistic for Sweden, each number appears in table (12) below represent 


degrees of freedom. 

-l) where X\-a ) is the theoretical value of Chi- 
of freedom. 

OF PARETO DISTRIBUTION: 
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the number of peoples for each Sweden cities at 2015 according to the population estimates, we regard each number as a 
random variables follows a Pareto distribution, so we write a computer program in Matlab R2013b to calculate (a,c) and 
using Chi-square test to find the values Q given by equation (7), where n=62. 


Table 12: Represents the Number of People in Each Cities in Sweden in 2015* 


2172770 

205162 

161938 

126841 

94077 

72259 

62172 

47494 

42392 

980347 

182315 

153952 

125113 

92132 

71213 

59089 

46024 

39826 

631873 

180221 

151025 

123464 

88690 

65492 

58855 

45451 

38927 

277441 

177396 

136167 

109125 

87476 

64936 

58669 

45180 

38565 

268207 

174306 

134669 

109073 

85342 

63848 

57391 

44215 

38016 

248558 

170848 

131661 

103979 

80282 

63428 

55576 

43922 

37417 

206161 

165153 

127376 

94845 

79316 

62176 

47758 

43425 



*The website of the Central Bureau of Statistics(SCB) Sweden ( http://www.scb.se/en / ) 


i- The Modified Moments Method: 


For this method, the estimation of (a,c)are: 

a = 1.3031 c = 369541.308 

The value of the test statistic for this variables is Q=10.19 with 5 deg. of freedom. This value of Q is compared 
with the tabulated critical values of a given sig. level of Chi-square test as shown in table (13) 


Table 13: Represents the Critical Values of a Given sig. Level of Chi-Square Test and Decisions for Sweden 


Chi-Squared 

Deg. of freedom 
Statistic 
P- Value 

5 

10.19 

0.07004 

Significance Level 

0.2 

0.1 

0.05 

0.02 

0.01 

Critical Value 

7.2893 

9.2364 

11.07 

13.388 

15.086 

Decision 

Reject 

Reject 

Accept 

Accept 

Accept 


ii- The Standard Bayes Method 

For this method, the estimation of (a,c)are: 

6 = 0.9984 6 = 37417 

The value of the test statistic for this variables is Q=14.279 with 5 deg. of freedom. This value of Q is compared 
with the tabulated critical values of a given sig. level of Chi-square test as shown in table (14) 

Table 14: Represents the Critical Values of a Given Sig. Level of Chi-Square Test and Decisions for Sweden 


Chi-Squared 

Deg. of freedom 
Statistic 
P- Value 

5 

14.279 

0.01393 

Significance Level 

0.2 

0.1 

0.05 

0.02 

0.01 

Critical Value 

7.2893 

9.2364 

11.07 

13.388 

15.086 

Decision 

Reject 

Reject 

Reject 

Reject 

Accept 
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The following chart shows the values of the p.d.f. at the values of random variables X l ,X~,,...,X 61 that are 
given in table (12). 
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Figure 9: The Values of p.d.f at Specific Values of Random Variables 

6-1-2 Maryland. 

For this application, each number in table (15)represent the number of peoples for each Maryland state and its 
Counties at 2015 according to the population estimators, we regard each number as a random variables follow a Pareto 
distribution, where n=24. 


Table 15: Represents the Number of People in Each Cities in Maryland State in 2015.* 


559600 

252000 

245600 

91650 

74650 

33900 

50150 

33250 

102950 

832050 

309050 

1036000 

157100 

30100 

103600 

39100 

26900 

52900 

168550 

625000 

900350 

113900 

151200 

20600 





*the source: united States Census Bureau (U.S. and World population Clock) http://www.census.gov/popclock/ 


i- The Modified Moments Method: 

For this method, the estimation of (a, C ) are: 


a = 1.0859 6 = 19810 

The value of the test statistic for this variables is Q=12.945 with 1 deg. of freedom. This value of Q is compared 
with the tabulated critical values of a given sig. level of Chi-square test as shown in table ( 1 6) 

Table 16: Represents the Critical Values of a Given Sig. Level of Chi-Square Test and Decisions for Maryland State 


Chi-Squared 

Deg. of freedom 
Statistic 
P- Value 

1 

12.945 

3.2083E-4 

Significance Level 

0.2 

0.1 

0.05 

0.02 

0.01 

Critical Value 

1.6424 

2.7055 

3.8415 

5.4119 

6.6349 

Decision 

Reject 

Reject 

Reject 

Reject 

Reject 
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iv- The Standard Bayes Method: 

For this method, the estimation of (a,c)are: 

a = 0.503 6 = 20600 

The value of the test statistic for this variables is Q=1.0217 with 2 deg. of freedom. This value of Q is compared 
with the tabulated critical values of a given sig. level of Chi-square test as shown in table (17) 


Table 17: Represents the Critical Values of a Given Sig. Level of Chi-Square Test and Decisions for Maryland State 


Chi-Squared 

Deg. of freedom 

Statistic 

P-Value 

2 

1.0217 

0.6 

Significance Level 

0.2 

0.1 

0.05 

0.02 

0.01 

Critical Value 

3.2189 

4.6052 

5.9915 

7.824 

9.2103 

Decision 

Accept 

Accept 

Accept 

Accept 

Accept 


The following chart shows the values of the p.d.f. at the values of random variables X l ,X 2 ,...,X 24 that are 
given in table (15). 
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Figure 10: The Values of p.d.f at Determined Values of Random Variables 


6-1-3 Alabama: 


For this application, each number in table (18) represent the number of peoples for each Alabama state and its 
Counties at 2015 according to the population estimators, we regard each number as a random variables follow a Pareto 
distribution, where n=67. 
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Table 18: Represents the Number of People in Each Cities and Counties in Alabama State * 


56223 

118324 

15354 

81996 

16832 

53171 

10629 

22487 

61337 

204543 

34637 

53269 

51229 

32157 

662177 

19246 

233033 

88886 

27182 

26756 

54720 

41463 

27284 

14421 

355475 

123102 

213585 

23367 

44236 

12886 

71996 

8722 

94572 

19819 

10084 

13411 

58466 

13426 

11232 

83426 

15807 

33234 

30549 

21327 

83262 

10711 

25765 

38505 

38281 

17760 

157737 

96108 

34108 

41656 

20347 

13817 

14196 

105019 

107711 

93509 

420180 

23185 

205652 

66163 

16510 

11154 

24354 







*the source: united States Census Bureau (U.S. and World population Clock) http://www.census.gov/popclock/ 


i- The Modified Moments Method: 

For this method, the estimation of (a , c ) are: 

a = 1.1324 6 = 8607 

The value of the test statistic for this variables is Q=47.662 with 5 deg. of freedom. This value of Q is compared 
with the tabulated critical values of a given sig. level of Chi-square test as shown in table (19) 


Table 19: Represents the Critical Values of a Given Sig. Level of Chi-Square Test and Decisions for Alabama State 


Chi-Squared 

Deg. of freedom 
Statistic 
P- Value 

5 

47.662 

4.1638E-9 

Significance Level 

0.2 

0.1 

0.05 

0.02 

0.01 

Critical Value 

7.2893 

9.2364 

11.07 

13.388 

15.086 

Decision 

Reject 

Reject 

Reject 

Reject 

Reject 


iv- The Standard Bayes method: 

For this method, the estimation of (a,c)are: 

a =0.6237 6 = 8722 

The value of the test statistic for this variables is Q=13.751 with 5 deg. of freedom. This value of Q is compared 
with the tabulated critical values of a given sig. level of Chi-square test as shown in table (20) 

Table 20: Represents the Critical Values of a Given Sig. Level of Chi-Square Test and Decisions for Alabama State 


Chi-Squared 

Deg. of freedom 

Statistic 

P-Value 

5 

13.751 

0.01727 

Significance Level 

0.2 

0.1 

0.05 

0.02 

0.01 

Critical Value 

7.2893 

9.2364 

11.07 

13.388 

15.086 

Decision 

Reject 

Reject 

Reject 

Reject 

Accept 


The following chart shows the values of the p.d.f. at the values of random variables X { ,X 2 , ...,X 61 that are 
given in table (18). 
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Figure 11: The Values of p.d.f at Determined Values of Random Variables 
6-2 Application (2): Tropical Cyclones: 

In this application we consider the tropical cyclones reaching TS (tropical storms) intensity or higher in Japan in 
2014 where each random variable is tropical cyclones with assuming that these variables have a Pareto distribution, so we 
write a computer program in Matlab R2013b to calculate (a,c) and using Chi-square test to find the values Q by equation 
(7) and The comparison between them. 

For this application the r.v. is the list of tropical cyclones reaching TS intensity or higher in 2014 in Japan as 
shown in table (21), where n=22. 


Table 21: Represents the Number Tropical Cyclones in Japan in 2014*. 


35 

35 

40 

70 

60 

45 

115 

45 

115 

35 

50 

100 

55 

75 

50 

110 

95 

40 

65 

40 

90 

110 







* the source: Annual Report on the Activities of the RSMC Tokyo - Typhoon Center 2014, Japan 
Meteorological Agency (2015) 


i- The Modified Moments Method: 

For this method, the estimation of ( a,C ) are: 

a = 2.0426 c = 34.221 

The value of the test statistic for this variables is Q=2.8081 with 2 deg. of freedom. This value of Q is compared 
with the tabulated critical values of a given sig. level of Chi-square test as shown in table (22) 


Table 22: Represents the Critical Values of a Given Sig. Level of Chi-Square 
Test and Decisions for Tropical Cyclones 


Chi-Squared 

Deg. of freedom 
Statistic 
P- Value 

2 

2.8081 

0.2456 

Significance Level 

0.2 

0.1 

0.05 

0.02 

0.01 

Critical Value 

3.2189 

4.6052 

5.9915 

7.824 

9.2103 

Decision 

Accept 

Accept 

Accept 

Accept 

Accept 
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ii- The Standard Bayes Method: 

For this method, the estimation of (a,c)are: 

a = 1.6202 c = 35 

The value of the test statistic for this variables is Q=0. 49268 with 2 deg. of freedom. This value of Q is compared 
with the tabulated critical values of a given sig. level of Chi-square test as shown in table (23) 


Table 23: Represents the Critical Values of a Given Sig. Level of Chi-Square 
Test and Decisions for Tropical Cyclones 


Chi-Squared 

Deg. of freedom 
Statistic 
P- Value 

2 

0.49268 

0.78165 

Significance Level 

0.2 

0.1 

0.05 

0.02 

0.01 

Critical Value 

3.2189 

4.6052 

5.9915 

7.824 

9.2103 

Decision 

Accept 

Accept 

Accept 

Accept 

Accept 


The following chart shows the values of the p.d.f. at the values of random variables X { ,X 2 X 22 that are 
given in table (21). 
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Figure 12: The Values of p.d.f at Specific Values of Random Variables 

7. CONCLUSIONS 
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From our theoretical and practical study, the following conclusions were pointed out: 

• B.M yield a small value of MSE comparing with M.M.M. 

• The bias of the estimators decrease as the sample size increase. 
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• The fit of tropical cyclones coincide with Pareto distribution for each level of significant, while, the fit of Sweden 
population coincide with Pareto distribution only for some levels of significant. 

• The goodness of fit test shows that the data of the population of Maryland state in U.S.A follows a Pareto 
distribution for each level of significant when the standard Bayes methods of estimation is used while the data 
does not appear to follow a Pareto distribution at any level of significant when the modified moment method of 
estimation is used. 

• The data of Alabama state follow a Pareto distribution only at 0.01 level of significant and the standard Bayes 
method is used. 

• The standard Bayes method of estimation is better than a modified moment in a sense of giving reasonable results 
for goodness of fit test. 

• The procedures studied in this paper can be extended to include the bivariate and generalized Pareto distribution. 
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